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III.  NEW  CONCEPTS  AND  EXPERIMENTS 


3.  1 Introduction 

This  volume  describes  some  new  hypotheses  and  experiments 
that  are  relevant  to  animal  echolocation  and  to  man  made  sonar  sys- 
tems. 

The  first  hypothesis  is  that  some  echolocating  animals  may 
use  very  wide  signal  bandwidths  in  order  to  compensate  for  a rela- 
tively small  array  that  has  only  two  transducers.  This  hypothesis 
depends  upon  the  assumption  that  a trade-off  exists  between  array 
size  and  signal  bandwidth.  The  existence  of  such  a trade-off  is 
investigated  by  devising  a new  ambiguity  function  that  applies  to 
range  and  angle  measurements,  i.  e.  , to  estimates  of  a target's 
position  in  space.  An  effective  beam  width  can  be  defined  in 
terms  of  the  new  ambiguity  function,  and  it  is  found  that  this  beam 
width  decreases  as  the  array  is  made  larger  (the  classical  result) 
and  also  as  signal  bandwidth  is  increased  (a  new  result).  A trade- 
off between  array  size  and  sigrial  bandwidth  therefore  does  exist  for 
a given  required  effective  beam  width.  This  result  supports  the 
hypothesis  that  the  use  of  large  signal  bandwidths  by  dolphins 
and  some  bats  may  serve  to  compensate  for  a limited  array  size. 

The  result  may  also  have  considerable  impact  upon  the  design  of 
future  man  made  sonar  systems. 

The  second  hypothesis  is  that  bats  (and  dolphins)  may  be 
able  to  estimate  cross -range  velocity  with  sufficient  accuracy  to 
compensate  for  evasive  maneuvers  by  moths  (and  fish).  This  hy- 
pothesis can  be  investigated  by  studying  the  dependence  of  expected 
error  in  a cross-range  velocity  estimate  upon  parameters  such  as 
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range,  interaural  distance,  and  signal -to -noise  ratio.  It  is  shown 
that  cross -range  velocity  estimation  is  indeed  a feasible  process 
for  targets  that  are  less  than  one  meter  from  the  animal,  where 
estimation  of  evasive  maneuvers  becomes  critical. 


The  third  hypothesis  involves  the  Johnson-Titlebaum  model 
of  echolocation  by  energy  spectrum  analysis.  The  hypothesis  is 
that  locally  optimum  detection  is  used  to  determine  the  presence 
or  absence  of  a target  by  investigation  of  the  energy  spectrum. 

The  energy  spectrum  is  obtained  by  Fourier  analysis  of  data  over 
a long  integration  time  that  includes  both  the  transmitted  signal 
and  the  target  echo.  The  locally  optimum  detector  correlates  sam- 
ples of  the  energy  spectrum  with  corresponding  samples  that  would 
be  measured  in  the  absence  of  noise.  The  noise-free  samples 
correspond  to  a ripple  or  sinusoidal  amplitude  modulation  of  the 
energy  spectrum,  where  the  period  of  the  ripple  is  determined 
by  target  range.  The  locally  optimum  test  for  a target  at  a se- 
quence of  hypothesized  ranges  therefore  involves  correlation  of 
. .........  • • 

the  energy  spectrum  with  a sequence  of  sinusoids,  i.  e.  , Fourier 

analysis  of  the  energy  spectrum.  This  hypothesis  suggests  that 
the  Johnson-Titlebaum  processor  is  particularly  sensitive  to 
background  noise  with  rippled  power  spectrum,  where  the  period  of 
the  ripple  corresponds  to  the  range  of  a target.  Since  a matched 
filter  is  comparatively  impervious  to  spectral  ripple,  we  have  the 
basis  for  a behavioral  experiment  that  can  distinguish  between  the 


Johnson-Titlebaum  and  matched  filter  models  of  animal  echo- 
location. 


Section  3.  2 describes  the  new  range-angle  ambiguity  func- 
tion and  the  bandwidth-array  size  trade-off.  Section  3.  3 presents 


a feasibility  argument  in  support  of  the  hypothesis  that  echolocating 
animals  can  measure  cross-range  velocity.  Section  3.4  discusses 
the  locally  optimum  detector  for  energy  spectrum  analysis  and  an 
experiment  to  test  the  Johnson- Titlebaum  echolocation  model. 

3.2  The  Trade-Off  Between  Array  Size  and  Signal  Bandwidth, 

and  Its  Possible  Utilization  in  Animal  Echolocation 

3.2.1  Brief  Summary  oi  Section  3.  2 

Target  position  estimation  in  radar  and  sonar  means  joint 
estimation  of  range  and  angle  in  the  presence  of  noise  and  clutter. 

The  global  behavior  of  a maximum  likelihood  position  estimator, 
and  the  clutter  suppression  capability  of  the  system,  can  be  written 
in  terms  of  a range-angle  ambiguity  function.  This  function  depends 
upon  signal  waveform  and  array  configuration,  i.  e.  , upon  both  tem- 
poral and  spatial  characteristics  of  the  system. 

Ambiguity  and  variance  bound  analysis  indicates  that  system 

• . . » • • • ' • « • 

bandwidth  can  often  be  traded  for  array  size,  and  direction- 
dependent  signals  can  be  used  to  obtain  better  angle  resolution  with- 
out increasing  the  size  of  the  array.  Wideband,  direction-dependent 
signals  (temporal  diversity)  can  be  traded  for  large  real  or  synthetic 
arrays  (spatial  diversity).  This  trade-off  is  apparently  exploited  by 
some  animal  echolocation  systems. 

The  above  insights  are  obtained  mostly  from  the  properties 
of  the  range -angle  ambiguity  function.  In  general,  an  appropriate 
ambiguity  function  should  be  very  useful  for  the  design  and  evaluation 
of  any  maximum  likelihood  parameter  estimator.  System  identification 
methods  and  radio  navigation  systems,  for  example,  can  be  optimized 
by  minimizing  the  volume  of  a multiparameter  ambiguity  function. 


3-2.2  Introduction 
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System  requirements  for  radar,  sonar,  or  diagnostic  ultra- 
sound include  parameters  such  as  array  size,  signal  bandwidth,  and 
processor  complexity  (measured  by  the  time-bandwidth  product  of 
the  signal,  data  storage  capacity,  and  processing  time).  System  per- 
formance involves  quantities  such  as  estimator  variance,  resolution, 
and  clatter  rejection  capability.  Performance  can  be  measured  in 
terms  of  Cramer-Rao  (CR)  bounds,  ambiguity  function  properties, 
and  signal-to-clutter  ratio.  If  these  performance  measures  can  be 
written  in  terms  of  array  size,  signal  bandwidth,  and  processor 
complexity  (system  requirements),  then  trade-offs  between  require- 
ments such  as  array  size  and  bandwidth  can  be  assessed. 

CR  bounds1,2  provide  a local  measure  of  estimator  accuracy 
for  particular  values  of  the  estimated  parameters.  In  terms  of  a 
hypothesis  test,  CR  bounds  are  useful  when  the  hypothesized  parame- 
ters are  nearly  equal  to  their  true  values,  i.  e.,  when  there  is  good 
prior  information  about  the  parameters  that  are  to  be  estimated.* 

CR  bounds  are  often  obtained  by  assuming  that  the  data  consists  of  a 
signal  that  depends  upon  the  unknown  parameters,  added  to 
white,  Gaussian  noise  (WGN).  If  the  data  actually  consists  of  addi- 
tional, spurious  signals  (clutter),  then  the  WGN  assumption  is 
invalid,  and  meaningful  CR  bounds  can  be  obtained  only  after  a more 
general  probability  distribution  has  been  worked  out. 

The  ambiguity  function3,4  illustrates  the  global  properties  of 
a maximum  likelihood  (ML)  estimator.  In  terms  of  a hypothesis 
test,  the  ambiguity  function  shows  the  response  of  the  estimator  to 
all  possible  values  of  the  hypothesized  parameters,  when  the  true 
parameters  are  specified.  Possible  confusion  between  different 
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parameter  values  can  thus  be  illustrated.  Like  CR  bounds,  the 
ambiguity  function  is  often  derived  under  a WGN  assumption,  but 
the  effect  of  spurious  echoes  can  be  analyzed  in  terms  of  this  idealized 
ambiguity  function,  and  its  utility  is  not  restricted  to  the  WGN  case.5  10 
In  view  of  these  properties,  it  is  not  surprising  that  ambiguity  analysis 
often  results  in  the  imposition  of  constraints  upon  a system  design  that 
minimizes  a CR  bound.  For  example,  the  ambiguity  function  gives  the 
accuracy  of  prior  knowledge  that  is  required  for  such  a system. 

We  shall  consider  properties  of  both  CR  bounds  and  the  am- 
biguity function  for  joint  range-angle  estimation.  Targets  that  are 
in  the  near  field  of  an  array,  as  well  as  the  far  field,  will  be  con- 
sidered, since  the  effective  size  of  synthetic  aperture  arrays  is  so 
large  that  the  far  field  assumption  is  often  violated.  Diagnostic 
ultrasound  arrays  for  examination  of  near- surface  effects  (breast 
and  parotid  tumors,  carotid  artery  wall  thickness)  will  also  some- 
times operate  in  a near-field  configuration. 

The  ambiguity  function  approach  to  radar /sonar  signal  design 
46 

is  well  established,  ^ and  this  paper  uses  a s^mjlar  method  for  com- 
bined signal  and  array  (temporal  and  spatial)  design.  Ambiguity 
analysis,  however,  has  apparently  been  neglected  for  other  parame- 
ter estimation  problems,  and  it  has  been  confined  to  the  radar /sonar 
literature.  One  purpose  of  this  paper  is  to  illustrate  the  advantages 
of  multiparameter  ambiguity  analysis45  in  a more  general  context. 

It  will  be  shown  that  the  volume  of  the  ambiguity  function  pro- 
vides a measure  of  the  extent  to  which  different  parameter  estimates 
can  be  separated  from  each  other  and  from  a background  that  may 
include  spurious  signals  as  well  as  noise.  By  far  the  best  known 
ambiguity  function  is  associated  with  delay  and  frequency  shift  mea- 
surements, using  narrowband  signals.3  The  volume  of  this  ambiguity 
function  depends  only  upon  signal  energy,  and  the  resulting  volume 
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invariance  for  energy-normalized  signals  can  be  interpreted  as  a 

radar  uncertainty  principle.27  The  volume  invariance  of  the  narrow- 

band  range -Doppler  ambiguity  function  is  more  the  exception  than 
J4  15 

the  rule,  however,  ’ and  the  analysis  in  this  paper  demonstrates 
that  minimization  of  ambiguity  volume  can  be  an  important  technique 
for  the  optimization  of  a parameter  estimator.  Minimization  of 
ambiguity  volume  should  be  especially  applicable  to  signal  design 
for  system  identification,  a branch  of  control  theory  that  has  much 
in  common  with  radar/sonar  design. 

Animal  sonar  or  echolocation  is  a biological  phenomenon 

that  is  especially  fascinating  to  radar /sonar  engineers.  One  of  the 

first  proposals  for  a man-made  sonar  system  was  suggested  in  1912, 

after  the  Titanic  disaster,  as  a simulation  of  "the  sixth  sense  of  the 

bat."26’28  Photographs  of  bats  in  pursuit  of  insects  among  foliage 

seem  to  indicate  a high  degree  of  clutter  suppression  capability,  and 

the  echolocation  performance  of  cetaceans  in  shallow,  reverberation- 

limited  environments  is  also  impressive.  Animal  sonars  operate  under 

a significant  constraint;,  their  receiving  "arrays"  consist  of  Qnly  two 

closely  spaced  elements.  Many  echolocating  animals,  however,  have 

evolved  extremely  wideband  systems,  and  the  highest  frequency  can 

■>o 

be  as  much  as  a factor  of  ten  larger  than  the  lowest  frequency.'  We 
know  that  increasing  a signal's  bandwidth  results  in  better  range  resolu 
tion3,4  and  better  target  discrimination  capability,30  55  and  that  large 
time -bandwidth  products  can  be  advantageous  for  wideband  velocity 
resolution36'37  or  for  Doppler  tolerance.38'40  We  also  suspect  that 
the  use  of  wideband  signals  may  be  associated  with  better  angle  mea- 
surements,11*21 and  that  echo-locating  animals  may  be  able  to  com- 
pensate for  limited  array  size  by  using  larger  bandwidths.  The 
following  analysis  indicates  that  this  suspicion  is  well-founded. 
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For  the  reader  who  is  interested  in  radar,  sonar,  diagnostic 
ultrasound,  and  animal  echolocation,  the  paper  should  provide  some 
useful  results  and  insights.  For  the  reader  who  is  interested  in 
other  parameter  estimation  problems,  the  specific  results  are  not 
so  important,  but  the  general  approach  illustrates  a useful  system 
design  philosophy. 

3.2.3  Problem  Formulation  and  General  Results 

In  this  section,  expressions  for  CR  bounds  and  the  ambiguity 
function  are  derived  for  the  position  estimation  problem.  The  results 
will  apply  to  all  planar  array  configurations.  No  assumptions  about 
array  shape  or  target  position  will  be  introduced,  except  for  the  simpli- 
fying assumption  that  the  target  lies  in  the  plane  of  the  receiving  array. 

Let  E(u)  be  the  Fourier  transform  of  an  echo,  measured  at 
the  target.  The  frequency  domain  response  of  the  k**1  array  ele- 
ment to  the  echo  is: 

Fk(u)  = \ E(tJ*  exPHWTk)  (1) 

where 

Ak  = complex  gain  of  the  kth  element  (2) 

= delay  from  target  to  kth  element  . (3). 

Let  L be  a straight  line  segment  between  the  elements  of  the  array 
that  are  furthest  apart.  The  range  parameter  t is  defined  as  the 
delay  between  the  target  and  the  center  point  P of  L,  i.  e.  , the 
point  that  lies  midway  between  the  extremities  of  the  array  (see 
Figure  1).  The  direction  parameter  0 is  defined  to  be  the  angle 
between  the  normal  to  L at  P and  the  line  segment  between  P and 
the  target.  Movement  of  the  target  counterclockwise  around  P 
causes  0 to  increase,  i.e.,  0 is  positive  in  the  counterclockwise 
direction. 
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The  position  of  each  array  element  can  be  defined  in  the 

same  way  as  the  target  position.  In  Figure  1,  the  k1*1  element 

position  is  described  in  terms  of  the  direction  parameter  0^  and 

the. distance  d from  P,  where  0 is  the  angle  between  the  nor- 
K K 

mal  to  L at  P and  the  line  segment  between  P and  the  array  ele- 
ment. From  Figure  1,  a general  expression  for  is 

Tk  = [t2  + (dk/c)2  - 2t  (dk/c)  <=°s(0k  - Q)]1^2  (4) 


and  depends  upon  the  unknown  parameters  t and  0 , as  well 

as  upon  the  known  parameters  d^  and  0^  . 

The  array  element  responses  in  the  absence  of  noise  can  be 

represented  by  a column  vector  m(T  , 0 ) of  complex  numbers. 

The  first  N elements  of  m(T,  9)  are  frequency  samples  of 

Fj(cu)  in  (1).  The  next  N elements  are  samples  of  F^(uo),  etc. 

In  representing  F,  (uu)  by  N frequency  samples,  it  is  implicitly 
^ th 

assumed  that  the  response  of  the  k element  is  time  limited. 
Although  a time  limited  signal  cannot  be  band  limited,  it  is  as- 
sumed that,  for  any  acceptable  mean  square  error  e>  0 between 
the  original  time  signal  and  the  time  signal  that  is  reconstructed 
from  N frequency  components,  there  exists  a bounded  number  N 
such  that  the  mean  square  reconstruction  error  is  less  than  e . 

For  an  integration  interval  T,  the  N discrete  frequency  samples 
are  separated  by  tt/T  rad/sec.  For  K array  elements,  m has 
KN  elements. 

In  Gaussian  noise,  a data  vector  £ is  observed,  where 
the  elements  of  r_  are  the  frequency  samples  of  Fj (uj),  F^ (uj),  . . . , 

F^Jiu)  when  noise  has  been  added  to  the  echo.  The  conditional 

^ H 

probability  density  function  of  the  noisy  data  is 
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p(r|T,0)  = (TTKN/r2|C|)_1  exp{-[r -m(T,e)]V  Q[r-m(T,0)]}  (5) 

where  C = Q * is  the  noise  covariance  matrix  and  the  asterisk  denotes 
a conjugate -transpose  operation. 

In  the  presence  of  multiple  targets  and  clutter  or  reverbera- 
tion, as  well  as  noise,  it  is  still  possible  to  use  a Gaussian  pdf  as  in 
(5),  provided  that  an  appropriate  covariance  matrix  is  determined.'2 
The  pdf  in  (5)  may  appear  to  be  unrealistic  for  sea  clutter,  since  the 
envelope  of  a Gaussian  clutter  process  is  Rayleigh  distributed,  while 
the  measured  pdf  of  envelope  detected  sea  clutter  echoes  is  often  more 
log-normal  than  Rayleigh. 41,42  The  empirical  log-normal  model  for 

the  clutter  envelope  indicates  that  the  underlying  pdf  of  the  time  domain 
clutter  process  decreases  more  slowly  in  the  "tails"  than  the  Gaussian 
distribution.  The  use  of  a frequency  domain  representation  in  (5), 
however,  implies  that  time  domain  echoes  are  effectively  passed 
through  a bank  of  narrowband  filters  with  impulse  responses  that  are 
T seconds  long.  The  sampled  responses  r_  of  these  filters  tend  to  be 
much  more  Gaussian  than  the  original  time  domain  distribution.43 

Maximum  likelihood  (ML)  estimates  of  range  and  angle  are 
the  values  of  t,0  that  maximize  p(r  Jt,  G)  in  (5),  for  a given  data 
vector  r.  The  asymptotic  variance  of  these  estimates  for  many 
observations  of  the  data  can  be  obtained  from  the  matrix  CR  bound.1’2 
The  CR  bound  is  defined  for  a particular  value  of  t,  0,  and  it  does 
not  take  account  of  the  possibility  that  two  very  different  parameter 
pairs  T1»®1  and  t2»  maY  Yield  nearly  equal  local  maxima  of 
p(r|T,  0).  The  ambiguity  function  is  a global  measure  of  ML 
estimator  performance  that  reveals  the  parameter  pairs  which  are 
most  likely  to  be  confused  with  one  another,  i.  e.  , which  are  asso- 
ciated with  an  ambiguous  interpretation  of  the  data.  We  shall 
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obtain  general  expressions  for  both  the  CR  bound  and  the  ambiguity 
function. 

When  range  T and  angle  0 are  both  unknown,  CR  bounds  for 
their  joint  estimates  are  determined  by  inverting  the  Fisher  informa- 
tion matrix,  J.  The  elements  of  J are 

Jn  = JTT  = -E[02/ar2)  lnp(r  |t,  0)] 

J12  = Jt0  = [(3^/3t  99 ) lnp(r|Y,  0)] 

(6) 

J21  “ JfiT  = -E[(3Z/30  9t)  in  p(r|  T,0)] 

J22  = Jee  = -E[(a2/902) ln  P(£.|  T.  0)]  • 

Substituting  (5)  into  (6),  we  obtain 

* 

J = 29m  /9t  Q 9m/9T 
TT  ~ - — 

J a = 2Re{9m  /9t  Q 9m/99} 

To  - 

J = 2Re{9m*/90  Q 9m/9x} 

0T  “ = ~ 

Jnn  = 29m  /90  Q 9m/90 

where  rn  = m (t,0)  are  the  expected  responses  in  (5),  i.  e.  , fre- 
quency samples  of  E^(w)  in  (1),  for  k = 1,  2,  . . . , K.  The  derivatives 
in  (7)  are  evaluated  at  the  particular  t and  0 values  that  correspond 
to  the  target's  true  position.  The  diagonal  terms  of  _J  * are  the  CR 
bounds  for  the  variance  of  the  range  and  angle  estimates  t and  0 , 


Var  (t-  t)  > [Jtt 
Var  ( 6 - © ) 2 [Jee 


-1 

(JtG  ^Gt  / ^00 

(JT0  J0T  / J n>]  • 


(8) 


In  order  to  obtain  simple  mathematical  expressions  for  the 
above  variances,  it  will  be  assumed  that 


Q 


(9) 


where  Nq  is  the  noise  power  spectral  density  and  1^  is  the  identity 
matrix.  This  assumption  means  that  the  noise  is  white  and  Gaussian, 
and  that  the  noise  at  each  of  the  K array  elements  is  statistically 
independent  of  the  noise  at  any  other  element.  The  CR  bounds  that 
are  obtained  by  using  (9)  are  generally  not  relevant  for  an  environment 
that  includes  multiple  targets  and  clutter  or  reverberation,  although  it 
may  be  possible  to  approximate  a diagonal  covariance  matrix  by  using  a 
large  integration  time  T for  the  computation  of  the  frequency  domain 
samples  £_.44  The  WGN  assumption  yields  idealized  bounds  that  con- 
vey only  qualitative  information  about  the  relation  between  system 
requirements  and  performance.  These  qualitative  insights,  however, 
can  be  very  important.  For  example,  CR  bounds  for  WGN  show  that 
the  accuracy  of  angle  estimates  are  not  completely  determined  by  the 
physical  beamwidth  of  the  receiver." 

Substituting  (9)  into  (7),  we  have 

/®  K 

u>2|E(w)|2  du.  J |Ak|2  OTk/0T)2  (10a) 
-co  k=l 
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(10b) 


JtG  = J0T 


2No_1  (T/ir)  / J |E(w)|  2 dw  £ (A^ 
^ k=l 


0-rk/9r)  Ork/D6) 


J60  = 2No 


/ K 

w2  |E(w)|  2 du>  ^ |Ak|  2 (0Tk/0  0)2  . (10c) 
-®  k=l 


In  (10),  t k is  given  by  (4)  and  is  shown  in  Figure  1.  Further  simpli- 
fications can  be  obtained  under  appropriate  assumptions  about  the 
relative  magnitudes  of  t,  t^,  and  dk/c.  Some  of  these  simplifica- 
tions will  be  discussed  in  the  sequel. 

A range-angle  ambiguity  function  is  obtained  by  studying  the 
behavior  of  the  conditional  probability  density  function  when  we  try  to 
determine  maximum  likelihood  estimates  by  trial-and-error.  This 
procedure  has  been  utilized  by  Urkowitz,  Hauer,  and  Koval'3  to  derive 
an  angular  ambiguity  function,  x(0»9tj)* 

The  ML  estimates  t,8  are  the  values  of  the  hypothesized  range 
and  angle  th>  9^  that  maximize  p[r  | t^,  0^]  or  that  minimize 

e = [r-m(TH,  0H)]*  Q [r  “ 5?(TH  * 0H>]  (HI 

where 

r_  = m (t,  0 ) + n 

and  where  n is  a vector  of  zero  mean,  complex  noise  samples. 

Assuming  that  (9)  is  true,  E{e}  is  the  mean  square  error  be- 
tween r and  m( t , 0 ).  Since 

H H 


X N 

m*  (t,  0 ) m(r,  0 ) = ^ |Fk(“n) 

k=l  n=l 


(12; 


we  see  from  (1)  that  1 1 rra  ( T , 0)||  does  not  depend  upon  T , and  it  is 

therefore  invariant  in  r and  0.  It  follows  that  E(e)  is  minimized 

when  x(T>  ttj*  0Tt)  *s  maximized,  where 

n H 

X(t,  th,  e,  Gh)  OC  Re  | m*  (r,  6)  Q m (th>  0h)} 

K _ 

H 5 ,Ak>  (1/2ir) 

I 

“fK-’Hk1]  d“ 

K 2 

= Z lAkl  R(Tk-THk'  <13> 

k=l 


where  R(t)  is  the  autocorrelation  function  of  the  echo  and,  from 

(4), 

THk  = K2  + ,dk/c)2  - 2 rH(dk/c)  COS  <V  eH)]i  • (U) 

The  trial-and-error  ML  estimator  can  be  conceptualized  as  a 

large  number  of  estimators,  each  matched  to  a different  set  of  values 

of  the  unknown  parameters.  The  ambiguity  function  x(t,  t , B,  8 ) 

H H 

is  the  output  of  these  estimators  when  a noise-free  signal  is  present  at 
the  input  to  the  system.  If  the  ambiguity  function  is  large  for  more 
than  one  set  of  hypothesized  parameter  values  (Tpj » f)^),  t^len  intro- 
duction of  noise  can  easily  lead  to  erroneous  estimates. 

The  range-angle  ambiguity  function  in  (13)  was  obtained  under 
the  WGN  assumption  (9).  Unlike  the  CR  bound  for  WGN,  the  ambi- 
guity function  in  (13)  can  be  applied  to  the  analysis  of  system  per- 
formance in  the  presence  of  clutter  or  reverberation  and  multiple 
5-10 

targets.  If  a point  target  is  at  position  t,6  and  a point  clutter 
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reflector  is  at  position  t , 9c  , the  output  power  of  a correlation  pro- 
cessor due  to  the  clutter  reflector  is 


clutter  response  (t,  0)  cc  I X(t,  t , 8,  0 )|  . 

i c c 


The  expected  output  power  due  to  a superposition  of  statistically 
independent  clutter  reflectors  is 

JT  oo 

Pc  =E  | clutter  response  (t,  0)jcc  j J p (T  , 0 ) 


-It  -■» 


|X(T,  T , 0,  0 )|  2 dT 
C C ' C 


where  p (tc»  0c)  is  a probability  density  function5  or  scattering 
function9  that  describes  the  distribution  of  clutter  in  range -angle 


space. 


The  correlator  response  to  the  target  is 


PT  cc  |X(t,t,0,0)| 


and  the  signal  to  clutter  ratio  is 


SCR  = P_/P_  - 
T C 


For  the  case  when  p(T  , 0 ) is  uniform,  e.  g. , when  there  is  no 

c c 

prior  knowledge  of  the  clutter  distribution,  we  have 


SCR  = P_/P„  = normalized  volume  of 
Cs  J. 


| x(t  , t^,  o,  e^j)  | 


One  of  the  well-known  properties  of  the  narrowband  range- 
velocity  ambiguity  function  is  volume  invariance,  i.  e. , the  same 


ambiguity  volume  is  obtained  for  any  energy  normalized  signal.  It 
would  appear  that  this  volume  invariance  property  does  not  apply 
to  other  ambiguity  functions,  however,  e.g.,  the  wideband  range- 
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velocity  ambiguity  function  ’ and  the  function  in  (13).  In  fact,  we 

shall  show  that  range -angle  ambiguity  volume  and  SCR  depend  upon 

the  bandwidth  of  the  radar /sonar  system. 

An  effective  beamwidth  can  also  be  written  in  terms  of  (16). 

For  a target  at  angle  0,  the  expected  response  from  clutter  at 

angle  0 , integrated  over  range,  is  the  bracketed  integral  in  (16). 
c 

If  this  integral  is  small  relative  to  P^,  then  the  effective  beam 
pattern  has  small  gain  at  0^  when  the  center  of  the  beam  is  aimed 
at  the  target.  For  uniformly  distributed  clutter. 

Effective  angular 
attenuation  = P 

T 

between  0 and  6 

c 


We  have  now  obtained  expressions  for  CR  bounds  and  the 
range-angle  ambiguity  function,  under  a WGN  assumption.  The 
general  target  and  array  description  in  Figure  1 has  led  to  compar- 
atively simple  notation  in  (10)  and  (13),  where  no  assumptions  about 
specific  geometries  have  been  introduced  (except  for  a planar  array 
with  a coplanar  target).  These  expressions  will  now  be  applied  to 


position  measurements  with  a linear  array. 

3.2.4  Linear  Array,  Far  Field 

For  a linear  array,  we  have  4>k  = ±"/2  for  all  K array  ele- 
ments. For  simplicity  of  notation,  we  shall  set  <t,  = -ir/2  and  we 

k 

shall  allow  d^  to  be  negative  as  well  as  positive.  In  this  case,  (4) 
becomes 
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(21) 


■’[ 


1 + t ^(d^/c)2  + 2t  * (d^/c)  sin  0 j 


If  the  target  is  in  the  far  field,  we  can  assume  that 


(d^/c)2  « t2 


and  (21)  becomes 


(22) 


= t + (d^/c)  sin  0 


(23) 


and 


i <vc)  i\|2  / i iv2- 

k=l  / k=l 


(28) 


Because  Jtq  = J0T,  the  product  Jtq  Jqt  in  (8)  is  always 
non-negative.  For  given  values  of  JTJ  and  Jq0,  and  for  Jtq  J0T  > 0, 
(8)  indicates  that  the  best  performance  is  obtained  when 


Jt0  J0t  = 0 - 


(29) 


According  to  (24c)  and  (28),  this  condition  is  obtained  in  WGN  when  the 
array  is  symmetric  about  its  midpoint,  i.  e.,  when 

d2  = "dr  d4  = "d3*  ’ • * 

(30) 

|a2!2  = (a^2,  Ia4I  2 = |a3I2,  . . . 


or 


d 


1 


= 0, 


IA31 2 = ia2i2,  ia5! 2 = |A4I2,  . 


(31) 


When  either  (30)  or  (31)  is  true,  (8)  becomes 


Var  ( r - x ) > [sNR  D^2  ] 
Var  (0  - 0)  s [sNR  dJ  D * 


(32) 

(33) 


Equality  in  (32)  and  (33)  is  asymptotically  obtained  for  a large  num- 
ber of  measurements,  if  a ML  estimate  is  used  in  WGN  (no  clutter), 
and  there  is  accurate  prior  knowledge  of  target  location. 
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What,  exactly,  is  meant  by  "accurate  prior  knowledge  of  target 

location"?  How  accurate  is  "accurate"?  Does  the  required  accuracy 

2 2 

of  prior  knowledge  change  when  D,.  and/or  D are  increased  in 

IV  A 

(32)  and  (33)?  These  questions  are  often  left  unanswered  in  the  deriva- 
tion of  parameter  estimation  methods  that  are  based  upon  CR  bounds, 
or  the  problem  is  side-stepped  by  assuming  that  the  estimates  are  con- 
sistent and  unbiased.  It  will  become  apparent,  however,  that  ambiguity 
analysis  provides  straightforward  answers  to  these  questions. 

The  bounds  in  (32)  and  (33)  can  be  minimized  by  using  two 
elements  at  opposite  ends  of  the  array  and  a narrowband  signal  with 
the  highest  allowable  frequency.  This  observation  follows  easily 
from  the  inequalities 


(34) 


(35) 


which  become  equalities  when  the  signal  energy  is  concentrated  at 
±W  and  only  two  array  elements  with  the  largest  d -value  are  used. 

K 

These  solutions,  however,  have  undesirable  ambiguity  properties,  as 
shown  below.  A ML  processor  should  be  evaluated  not  only  by  CR 
bounds  but  also  by  ambiguity  function  analysis. 

Although  significant  caveats  exist  for  quantitative  analysis 
of  (32)  and  (33),  the  bounds  are  qualitatively  significant  because  they 
illustrate  a basic  interdependence  between  signal  (temporal)  design 
and  array  (spatial)  design.  This  interdependence  will  become  even 
more  apparent  when  ambiguity  functions  are  considered. 


J 
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The  range-angle  ambiguity  function  for  targets  in  the  far  field 
of  a linear  array  is  obtained  by  substituting  (23)  into  (13), 

K 

X(t , th*  e»  ^ *Ak*2  R [T  ‘ th  + (dk/c)  (sin  6 - sin  0H)J 


X(t-t  sin  6 - sin  0 ) . 
H H 


If  0 and  8 are  less  than  30°, 
H 


X(T-  TH'  e'  °H)  ”2  |Ak|2R  [T-  TH  + (dk/c)(e-  V] 


= X(At,  A0 ) . 


When  the  angle  hypothesis  is  correct,  we  have  0 = 0 or 

A0  = 0,  and  (36)  is  proportional  to  R(At),  the  echo  autocorrelation 
function.  Additional  array  elements  do  not  affect  the  structure  of 
X(At,  0),  and  this  structure  includes  many  undesirable  sidelobes 
when  the  signal  is  narrowband.  Figure  2 shows  X(At,  A0)  for  a 
narrowband  pulse  and  with  two  hydrophones,  i.  e.  , for  the  "optimum" 
system  design  that  is  obtained  from  CR  bounds  as  in  (34)  and  (35). 

For  a single  target  in  a clutter-free  environment.  Figure  2 
indicates  that  an  accurate  position  estimate  can  be  obtained  only  if 
prior  knowledge  restricts  the  search  area  to  the  central  lobe  of  the 
ambiguity  function.  The  central  lobe  becomes  more  narrow  as  the 
signal  frequency  and  the  distance  between  hydrophones  increase,  so 
a more  accurate  estimate  implicitly  requires  better  prior  knowl- 
edge of  the  estimated  parameters.  Under  clutter-free  conditions, 
one  might  then  expect  the  two-hydrophone,  high  frequency  system  to 


t.W. 


TWO  HYDROPHONES,  NARROWBAND  PULSE 


I 


; j 

j 

j 

i 

1 
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emerge  as  the  solution  of  a sequential  beam  forming  procedure.16 
Ironically,  such  a final  result  would  be  a poor  starting  point  for  an 
adaptive  beam  former,  if  the  system  begins  without  prior  knowledge 
of  t and  0.  The  optimality  of  the  two  hydrophone  configuration 
depends  upon  accumulated  information  (prior  probability  densities 
with  small  variance),  and  the  result  is  useless  without  this 
information. 

Far  field  range  ambiguities  depend  upon  the  behavior  of 
R(At),  the  echo  autocorrelation  function.  R(At)  becomes  more 
impulse-like  when  the  bandwidth  of  the  echo  is  increased,  and  range 
ambiguities  can  be  reduced  by  using  a wideband  signal.  In  the  two 
hydrophone  case,  the  effect  of  switching  to  a wideband  signal  is  illus- 
trated in  Figure  3. 

Figure  3 illustrates  a decomposition  of  the  ambiguity  function 
into  a sequence  of  shifted  autocorrelation  functions  when  A0  becomes 
large.  This  effect  is  predicted  by  (37).  The  width  of  each  auto- 
correlation function  is  approximately  1/B  sec,  where  B is  the  echo 
bandwidth.  If  (d^/c)  A0  is  greater  than  1/B  for  all  k in  (37),  then 
a constant-A0  profile  of  the  ambiguity  function  is  a sequence  of  K 

2 

nonoverlapping  autocorrelation  functions,  each  with  amplitude  lA^l  . 
Decomposition  occurs  when  |A0|  > A0q,  where 

A0_  = (c/B)  / min  (d  ) 

0 k K 

(38) 

spatial  width  of  echo  autocorrelation  function 
minimum  distance  between  hydrophones 

The  ratio  in  (38)  can  be  small  if  very  large  bandwidths  are  used. 

In  the  case  of  sonar,  if  B = 50  kHz  and  array  elements  are  1 meter 
apart,  then  A0^  = .03  radian  (less  than  2°). 
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For  | A8l  > A0Q,  it  is  easy  to  obtain  some  fundamental 
properties  of  the  ambiguity  function.  Since  X(At,  A0)  for  A0  constant 
and  |A0|  > A0q  is  a sum  of  nonoverlapping  autocorrelation  functions, 
we  have 


max  X(At,  A0)  / X(0,  0) 
At 


|A0|  > A0 


K 

= max  |Ak|2R(0)/j  |Ak|2R(0) 

k / i._i 


> 1/K 


with  equality  when  |Ak(  is  the  same  for  all  k-values  (uniform 
weighting).  The  sidelobe  level  for  | A0|  > A0^  can  be  reduced  only 
by  using  a larger  array.  For  example.  Figure  4 shows  X(At,  A0) 
for  K = 5,  rather  than  for  K = 2,  as  in  Figure  3. 

In  terms  of  signal-to-clutter  ratio,  (39)  describes  the  maxi- 
mum response  of  the  system  to  a single  clutter  reflector  that  is 
A0  radians  from  the  target,  divided  by  the  response  to  the  target 
itself.  Eq.  (39)  is  thus  a measure  of  worst-case  SCR,  where  the 
clutter  is  restricted  to  a single  point  in  t,8  space.  For  uniformly 
distributed  clutter,  SCR  for  a given  value  of  A0  is  related  to  the 
effective  beamwidth  of  the  system,  as  described  by  (20).  Eq.  (20) 
gives  the  effective  angular  attenuation  between  the  center  of  the  beam 
and  uniformly  distributed  clutter  that  is  A0  radians  from  the  center. 
For  |a©|  > A0Q,  (20)  becomes 


WIDEBAND  SIGNAL,  FIVE  UNIFORMLY 
SPACED  HYDROPHONES 


J 


||a0|  > A60  ' 


K 


K 


- i Kl4  f |R(T)|2dT  / |£  |\|2K10>' 

k=l  -®  / k=1 


By  the  Cauchy  and  Schwarz  inequalities. 


and 

Effective  attenuation  (A0) 


> 1 /KB  , 


i |2 

with  equality  when  jA^)  = constant  (uniform  weighting  in  space) 
and  when  |E(u>)|2  = constant  over  B (uniform  weighting  in  frequency). 
These  equality  conditions  are  the  converse  of  the  equality  conditions 
for  (34)  and  (35),  which  minimize  CR  bounds  by  using  the  ultimate  in 
non-uniform  weighting. 

For  uniformly  distributed  clutter,  the  maximum  sidelobe 
level  in  (39)  is  not  as  relevant  as  the  effective  attenuation  measure 
in  (20)  and  (41).  In  order  to  produce  a small  effective  beamwidth 
for  uniformly  distributed  clutter,  the  product  KB  should  be  large, 
and  system  performance  depends  equally  upon  bandwidth  and  array 
size.  For  a given  array  size,  clutter  suppression  performance  can 
t be  improved  by  increasing  the  bandwidth  of  the  system. 

A more  general  measure  of  signal-to-clutter  ratio  is  the 
normalized  ambiguity  volume  in  (18)  and  (19).  It  is  shown  in  the 
Appendix  that,  when  |E(u>)|  is  constant  over  a bandwidth  B, 

SCR  > B/(2tt)2  . (42) 

For  uniformly  distributed  clutter,  SCR  increases  with  bandwidth. 

An  upper  bound  can  be  obtained  for  a two  element  system 
which  uses  a maximum  wavelength  \ and  a bandwidth  B such  that 

X S the  distance  between  receiving  elements 
B>  one  octave. 

These  conditions  hold  for  many  animal  sonar  systems,  e.  g.  , the 
bottlenosed  dolphin  (Tursiops  truncatus)  and  the  large  brown  bat 
(Eptesicus  fuscus).  For  these  animals, 

SCR  < 2B  (43) 
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when  clutter  is  uniformly  distributed.  Many  animal  echolocation 
systems  have  apparently  compensated  for  their  restricted  array 
size  by  using  large  bandwidths. 


In  summary,  analysis  of  the  linear  array,  far  field  case  has 
demonstrated  that  optimization  on  the  basis  of  CR  bounds  alone  is 
not  a desirable  procedure  unless  the  environment  is  clutter-free  and 
there  is  good  prior  knowledge  of  the  estimated  parameters.  The 
required  accuracy  of  the  prior  knowledge  can  be  deduced  from  an 
appropriate  ambiguity  function.  Qualitative  interpretation  of  CR 
bounds  yields  important  insight  into  the  fundamental  trade-offs 
between  signal  and  array  parameters.  These  trade-offs  can  be 
further  investigated  by  computing  bounds  on  signal-to-clutter  ratio, 
using  the  range-angle  ambiguity  function.  For  uniformly  distributed 
clutter,  these  bounds  depend  upon  bandwidth  alone  or  upon  the  pro- 
duct of  bandwidth  and  array  size.  Array  size,  independent  of  band- 
width, becomes  important  when  the  clutter  consists  of  a single 
reflector  with  particular  values  of  t and  6 . 


3.2.5  Linear  Array,  Near  Field 

For  synthetic  aperture  systems  and  for  some  diagnostic 
ultrasound  applications,  targets  are  often  in  the  near  field  of  the 
array.  The  elements  of  the  Fisher  information  matrix  are  in  this 
case  obtained  by  substituting  (27)  into  (10).  The  results  are 


JTT  = SNR 


D-2 1 •„*  Ki2  / 1 h 

k=l  / k=l 


(44a) 


K / K 

J6e  = SNRD„2co92e  £ Pk2  <VC)2K|2  / £ KI2  <44b» 

W = 1 / lc  = 1 


/ 


B = t y t = to  center  of  array 

Pk  k , . t , th  , * (4bb> 

delay  to  k element 

The  expression  for  J TT  indicates  that  the  array  configuration  is 
important  for  near-field  range  estimation,  a condition  that  did  not 
exist  for  far  field  measurements.  As  in  the  far  field  case,  J00 
depends  upon  both  temporal  and  spatial  parameters.  If  the  gains 
Ak  are  to  be  adjusted  so  as  to  maximize  JTT  and  J00,  the  optimum 
weights  will  depend  upon  target  position.  This  observation  suggests 
that  the  gains  A^  are  no  longer  decoupled  from  the  position  hypothesis 
t„i  and  the  gains  should  change  as  the  environment  is  scanned. 

From  (45),  the  element  that  is  closest  to  the  hypothesized  target 
position  should  have  the  largest  gain.  The  expression  for  JT0  J0t 
indicates  that  the  elimination  of  range-angle  error  coupling  is  more 
difficult  to-accomplish  in  the  near  field  case.  This  difficulty  arises 
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because  differences  in  the  near-field  delays  depend  upon  range 
as  well  as  angle. 

The  expression  for  the  near-field  ambiguity  function  can  be 
simplified  if  in  (27)  can  be  written 

« r|l  + (1/2)  [(dk/cT)2  + 2 (dj^/cT)  sin  e]  } . (46) 

If  the  term  in  square  brackets  is  less  than  unity,  the  error  in  the 
above  approximation  is  less  than  6.4%.  We  are  thus  assuming  that 
the  distance  from  the  target  to  the  center  of  the  array  is  larger  than 
the  array  itself.  Substituting  (46)  into  (13), 


x<T'Ve>v  = 2 Ki2  rKat  +(dk/c) 


(sin  0 - sin 

where 

sk  = 1 - [(dk/c)2  / (2tth)]  • (48) 

The  sum  in  (47)  can  have  an  effect  that  is  similar  to  broadening 
the  bandwidth.  When  0 = 0H>  the  near-field  ambiguity  function  depends 
upon  a superposition  of  weighted  autocorrelation  functions,  where  each 
autocorrelation  function  is  scaled  by  a factor  s^.  If  the  echo  autocor- 
relation function  has  sidelobes  (local  maxima  at  AT  t 0),  the  sum  of 
scaled  functions  can  have  a smaller  sidelobe  level  than  a sum  of  unsealed 
functions,  since  the  maxima  for  the  scaled  functions  at  AT  / 0 occur 
at  different  locations  for  different  scale  factors,  s^  . We  therefore 
suspect  that  the  range  resolution  capabilities  of  near  field  systems 
(e.  g.  , synthetic  aperture)  will  often  be  better  than  we  would  predict 
from  the  signal  autocorrelation  function. 


When  the  distance  from  the  target  to  the  array  is  smaller  than 

the  size  of  the  array  itself,  the  best  way  to  obtain  a picture  of  the 

range-angle  ambiguity  function  is  probably  to  evaluate  (13)  by  means 

of  a computer.  A function  of  t , 0 can  be  displayed  for  any  parti- 

cular  values  of  t and  0.  The  effect  of  sensor  positions,  gains,  and 

signal  bandwidth  upon  ambiguity  volume  can  then  be  empirically 

determined.  If  ambiguity  volume  is  used  as  a measure  of  system 

performance  (SCR),  then  the  optimum  signal/array  design  can  be 
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obtained  by  gradient  techniques. 


3.2.6  Direction- Dependent  Signals  and  Array  Elements 

The  effect  of  transmitting  an  angle-dependent  signal  is  to 

replace  |e(w)|^  in  (13)  with  E(w,  0)  E^(w,  0 ).  If  the  inverse  Fourier 

H 

transform  of  E(w,  0)  E (w,  0 ) has  small  maximum  amplitude  for 

H 

0 0,  the  ambiguity  function  will  be  rapidly  attenuated  for  A0  ^ 0. 

H 

The  most  common  method  of  achieving  direction  dependence 

is  through  the  use  of  a narrow  physical  beamwidth.  Other  methods 

include  the  use  of  a frequency-steered  array,  a dispersive  lens  (for 

sonar),  and  movement  of  the  array  with  respect  to  the  environment, 

so  that  an  angle-dependent  Doppler  history  is  obtained  from  each 

target.  The  latter  method  is  employed  in  synthetic  aperture 
a. 17,18  i _ j j , l i j ~ i l j 


systems . 


Angle-dependent  Doppler  shifts  could  also  be  used  in 


ultrasonic  blood  flow  measurements,  where  movement  of  the  target 
would  replace  movement  of  the  array. 

Direction-dependent  array  elements  can  also  increase  angle 
resolution.11  To  include  the  effect  of  a direction  dependent  transfer 
function,  the  gain  A^  in  (2)  is  replaced  by  A^(0).  More  generally, 
we  can  take  account  of  the  temporal  impulse  response  of  the  element 
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by  using  A^(0,  w).  The  quantity  |a^|  in  (13)  is  then  replaced  by 

A,  (0,  w)  A.  (0„,  lj) • The  effect  is  again  to  reduce  the  level  of 
k k H 

X(At,  A0)  for  nonzero  A0. 

A direction-dependent  receiving  element  that  makes  use  of 

19 

multiple  reflections,  as  in  the  human  pinna,  imposes  a direction- 

dependent  convolutional  code  upon  received  wideband  waveforms.  If 

the  receiver  can  decode  the  resulting  signal,  good  angle  resolution 

can  be  obtained  with  a very  compact  physical  array.  The  trade-off 

here  is  between  array  size  and  processor  complexity. 

It  is  likely  that  animal  echolocation  systems  exploit  direction 

dependence  of  transmitted  and  received  signals."  A transmitted 

dolphin  echolocation  pulse  has  different  structure  when  it  is  observed 
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at  different  angles  relative  to  the  animal,  and  the  external  ears 
of  many  bats  are  capable  of  imposing  a direction-dependent  code 
upon  received  signals. 

3.2.7  Extended  Targets  and  System  Identification 

The  formulation  in  (11)  to  (13)  can  produce  a generalized 
ambiguity  function*3  for  any  ML  estimation  problem,  and  it  would 
seem  that  the  ambiguity  function  should  not  be  restricted  to  radar/ 
sonar  applications.  Ambiguity  analysis  should  be  especially  bene- 
ficial for  system  identification,  where  CR  analysis  has  already  been 
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applied  to  the  design  of  probing  signals.  The  shortcomings  of 
signal  derivations  that  are  based  upon  CR  bounds  alone  have  already 
been  discussed.  It  remains  to  demonstrate  the  utility  of  ambiguity 
analysis  for  system  identification.  The  advantages  of  ambiguity  anal- 
ysis should  be  obvious  if  we  can  find  a system  identification  problem 
that  involves  radar/sonar  measurements.  Parameterization  of 
extended  targets  is  ideal  for  this  purpose. 


I 

An  extended  target  can  sometimes  be  described  as  a distri- 
bution of  point  reflectors  or  highlights  in  range-angle  space.  Most 
diagnostic  ultrasound  and  synthetic  aperture  imaging  systems  are 
based  upon  such  a description.  In  estimating  the  position  of  each 
highlight,  reflections  from  the  other  highlights  can  be  regarded  as 

clutter.  We  have  seen  that  CR  bounds,  derived  under  WGN  condi-  j 

tions,  are  not  applicable  to  this  situation.  An  impulse-like  ambiguity 

function,  however,  will  be  capable  of  resolving  the  target  into  its 

separate  highlights,  or  at  least  of  determining  the  target  reflectivity 

within  a small  range-angle  cell. 

An  extended  target  can  also  be  characterized  as  a distributed 
parameter  system.  The  impulse  response  is  a function  of  time  (range) 
for  a given  angle,  and  the  impulse  response  changes  with  angle.  The 
system  is  to  be  parameterized  by  the  locations  (i.  e.  , the  t,0  values) 
of  large  local  maxima  in  the  impulse  response  (i.  e.  , the  highlights). 

Each  of  these  position  parameters  is  to  be  estimated  in  the  presence 
of  noise  and  spurious  signal  components  (clutter)  that  are  part  of  the 
system's  response  to  the  probing  signal. 

The  system  identification  problem  for  an  extended  target  is 
thus  equivalent  to  the  estimation  of  a sequence  of  t,0  values  in  clut- 
ter. The  ambiguity  function  is  a valuable  aid  in  the  design  of  signals 
and  filters  that  can  perform  the  required  identification.  Exclusive 
use  of  CR  bounds  for  design  of  the  probing  signal  for  this  identification 
problem  can  result  in  an  ambiguity  function  with  large  volume  or  large 
sidelobes,  and  undesirable  interaction  between  parameter  estimates  — — i 

will  result.  The  concept  of  resolution  in  radar /sonar  performance 
translates  into  parameter  separability  or  increased  observability 
in  control  theory.  In  a system  identification  context,  the  volume  of 
the  ambiguity  function  is  a measure  of  the  extent  to  which  each  parame- 
ter estimate  is  unaffected  by  the  remaining  system  parameters. 
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3.2.8  Conclusion 

A range-angle  ambiguity  function  can  be  used  to  predict 
signal -to -interference  ratios  for  high-resolution  radar/sonar  sys- 
tems. The  expected  clutter  response  will  often  depend  upon  the 
volume  under  the  sidelobes  of  the  ambiguity  function,  rather  than 
the  height  of  these  sidelobes.  For  a given  array  size,  the  sidelobes 
are  made  "thinner"  by  using  a wideband  signal,  the  volume  is  thus 
reduced,  and  signal-to-interference  ratio  is  increased. 

Direction-dependent,  wideband  echoes  can  result  in  an  impulse- 
like range-angle  ambiguity  function.  In  synthetic  aperture  systems, 
the  direction  dependence  is  obtained  by  moving  the  radar  or  sonar 
relative  to  the  environment.  Other  methods  that  do  not  depend  upon 
movement  can  perhaps  be  used  to  increase  the  rate  at  which  a given 
area  can  be  mapped.  For  example,  one  can  use  an  angle-dependent 
pulse  with  large  time -bandwidth  product,  and  this  signal  could  be 
transmitted  in  all  directions  simultaneously. 

It  -would  seem  that  the  range-angle  ambiguity  function  is  a 
radar /sonar  counterpart  of  the  point  spread  function  that  is  used  to 
define  visual  acuity.  It  is  therefore  a useful  concept  for  sonar  sys- 
tems that  are  attempting  to  "see  with  sound.  " The  trade-off  between 
bandwidth  and  array  size,  which  has  been  obtained  from  properties 
of  the  range -angle  ambiguity  function,  would  seem  to  be  important 
for  echolocating  animals  that  use  wideband  signals. 

Cramer-Rao  bounds  are  helpful  indicators  of  qualitative  inter- 
dependencies between  temporal  and  spatial  parameters.  The  exclu- 
sive use  of  CR  bounds  for  signal  and  array  synthesis,  however,  can 
lead  to  results  which  are  only  optimum  in  the  immediate  neighbor- 
hood of  a particular  position,  i.  e.  , with  accurate  prior  knowledge 
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of  t and  0 . The  global  behavior  of  the  maximum  likelihood  estimate 
is  portrayed  by  the  ambiguity  function,  which  provides  a more  reliable 
functional  for  synthesis  of  signals  and  array  configurations.  This 
observation  can  be  generalized  to  any  parameter  estimation  problem, 
and  it  appears  to  be  especially  relevant  to  the  design  of  probing  signals 
for  system  identification. 

Multiparameter  space-time  ambiguity  analysis  can  be  used  to 
determine  the  best  array  locations  for  passive  sonar  systems,  and  the 
approach  is  also  applicable  to  signal  design  and  transmitter  place- 
ment for  radio  navigation. 

In  many  cases,  ambiguity  volume  depends  upon  the  properties 
of  the  estimation  device,  and  a small  volume  is  indicative  of  parameter 
separability.  Computer  minimization  of  the  volume  of  a multi- 
dimensional ambiguity  function  should  result  in  an  optimum  estimating 
device,  e.  g.  , an  optimum  probing  signal  for  system  identification. 
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3.  2.  10  Appendix  to  Section  3.  2 

Signal  to  clutter  ratio  is  defined  as  the  ratio  P _/P  , where 

T C 


PT  = IX(t,  t,  0,  6)1 


(A  - 1 ) 


pc  “ 


oo  TT 

J drc  j dec  p(tc,  ec)  ix(t,  Tc,e,  0C)|2  . 


(A- 2) 
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In  (A-2),  P(TC»  9c)  is  the  clutter  probability  density  function  or  a 
normalized  version  of  the  clutter  scattering  function,  and 

K * 

X(t,tc,6,0c)  = ( 1 / 2 it)  ^ lA^I  / IE(w)|^  exp|jwj^(d^/c) 

k=l  -oo 

(sin  0 - sin  0 ) + t - t J j du 


X(  t - t , sin  0 - sin  0 ) 
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(A- 3) 


for  a far -field  condition. 

For  uniformly  distributed  clutter. 
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Performing  the  t integration,  we  have 
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(1/2tt)  J exp  [-j(coi  - w2)  J 


d t = 6(w_  - w,  ) 
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and 


The  imaginary  part  of  the  0 -integral  is  zero,  i.e., 

c 

IT 

— (1/2tt)  J sin  [(<j/c)  (dk  - d^)  sin  ©c  ] d0c  2 0 


because  the  integrand  is  an  odd  function  of  0^  and  the  limits  of  inte- 
gration are  symmetric  about  0 =0.  The  real  part  of  the  0 inte- 
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gral  is 
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where  Jq(‘)  is  a Bessel  function  of  order  zero.  It  follows  that 


= 2 |Ak!2  lA"|Z  /lE(“>|4  Jot(“/c>  (dk-d„»] 

k,  n 'Lx> 


exp  j (co/c)  (dk-dn)  sinej  dco  . 
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The  expected  clutter  response  for  uniformly  distributed  clutter  is 
therefore  a function  of  0,  the  angle  between  a line  drawn  from  the 
center  of  the  array  to  the  target  and  a line  that  is  normal  to  the 
array.  An  upper  bound  on  P can  be  obtained  by  noting  that 


J0  [Vvc)  (dk-  dn)]  exp  [j  (co/c)(dk- dn)  sine] 


S 1 , 


and  therefore 
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From  (A-l)  and  (A-3), 
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(A -8) 
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If  |E(to)|  is  constant  over  a bandwidth  B,  then,  from  (40), 

PT/PC  > B / (2tt)2  . 

For  uniformly  distributed  clutter,  the  lower  bound  (A- 10)  on  signal- 
to-clutter  ratio  will  increase  as  bandwidth  increases. 

An  upper  bound  for  SCR  can  be  obtained  when  (i)  the  maximum 
echo  wavelength  is  less  than  the  minimum  distance  between  hydro- 
phones, and  (ii)  the  echo  spectral  magnitude  is  smooth  and  has  at 
least  an  octave  bandwidth.  The  first  condition  means  that2S 

i 

Jq  (x)  ~ (2/ttx)2  cos  (x  - 7r/4) 

where 

x = (w/c)  (d  - d ) . 

k n 

The  error  in  the  approximation  is  less  than  3.2 % for 
X > 2 7T 

or 


> X, 

which  is  condition  (i). 

The  second  condition  implies  that,  for  k ^ n. 
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because  the  smooth  non-negative  function  |e(co)|^  is  multiplied  by 
a periodic  function  with  at  least  one  full  oscillation.  Substituting 
(A- 11)  into  (A- 6), 


P /Pi 
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and  from  (40) 
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The  approximation  in  (A-ll)  must  be  used  carefully,  since 
the  sum  of  the  left-hand  side  over  a large  number  of  elements  may 
not  be  small  relative  to  the  right-hand  side,  even  though  each  term 
of  the  sum  is  small.  This  reservation  is  inconsequential  for  studies 
of  animal  echolocation,  where  K = 2.  The  two  conditions  that  lead 
to  (A-13)  are  satisfied  by  many  wide-band  animal  echolocation  systems. 
For  the  Atlantic  bottlenose  dolphin,  Tursiops  truncatus,  echolocation 
pulses  cover  a frequency  range  between  about  20  and  150  kHz.20  The 
maximum  wavelength  in  water  is  7.  5 cm,  which  is  shorter  than  the 
distance  between  the  ears.  The  large  brown  bat,  Eptesicus  fuscus, 
has  a cruising  pulse  that  covers  the  frequency  range  between  about 
25  and  65  kHz.26  The  maximum  wavelength  in  air  is  1 . 3 cm,  which 
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is  again  shorter  than  the  distance  between  the  ears.  Although  many 
animals  are  restricted  to  two  array  elements,  they  can  obtain  good 
SCR  by  using  very  wide  bandwidths. 
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3.  3 Measurement  of  Tangential  Velocity  and  Evasive  Maneuvers 
3.3.1  Introduction  to  Section  3,3 

The  utilization  of  evasive  maneuvers  by  moths  is  well  docu- 
mented for  bat  echolocation Noctuid  moths  that  are  10  to  100  feet 
away  from  a bat  generally  fly  away  from  the  sound  source,  and  those 
that  are  close  enough  to  be  tracked  display  erratic  flight  patterns. 
Tangential  motion  (i.e.,  motion  perpendic  lar  to  the  direction  of 

propagation)  is  an  optimum  maneuver  for  the  differential  game 

2 

known  as  the  homicidal  chauffeur  problem. 

Although  evasive  action  by  fish  is  less  well  documented,  the 
use  of  a sudden  tangential  motion  to  escape  a threat  at  close  range  is 
familiar  to  any  skin  diver  or  aquarium  hobbyist.  If  a dolphin  is  to 
catch  prey  at  night  or  in  very  turbid  water,  high  survival  value 
would  be  attached  to  the  detection  of  tangential  motion,  i,e.,  motion 
that  is  normal  to  the  dolphin's  own  velocity  vector,  at  close  range. 

What  is  the  best  way  to  measure  tangential  velocity,  and  how 

might  it  be  measured  in  animal  sensory  systems?  It  would  appear 

that  a maximum  likelihood  estimate  of  cross-range  velocity  can  be 

obtained  by  frequency  shifting  (or  time  scaling)  the  input  to  one 

ear  and  correlating  this  shifted  (or  scaled)  signal  with  the  input  to 

the  other  ear.  Maximum  correlator  response  will  occur  when  the 

Doppler  shift  has  been  correctly  hypothesized.  This  frequency  shift 

and  cross  correlation  is  analogous  to  the  time  shift  and  cross  corre- 

lation  that  has  been  suggested  by  Licklider  and  Sayers  and  Cherry 

as  a model  for  binaural  interaction.  An  equivalent  operation  can  also 

be  implemented  by  the  equalization-cancellation  model  that  has  been 

6 

proposed  by  Durlach. 


If  we  were  to  simultaneously  estimate  interaural  delay  (for 
angle  measurement)  and  interaural  Doppler  shift  (for  tangential 
velocity  measurement)  then  the  output  u^(t)  of  one  ear  should  be 
frequency  shifted  and  delayed,  and  then  correlated  with  the  output 
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u^(t)  of  the  second  ear.  The  correlator  response  is  a function  of 
the  frequency  shift,  0,  and  the  delay,  T.  This  response  is  the  cros 
ambiguity  function  of  the  echo,  Xu  Uo(t,  0). 

* ” X c* 

There  are  some  interesting  theoretical  implications  to  the 
conjecture  that  the  cross-ambiguity  function  xUjU2(t,  p)  is  formed 
by  the  auditory  system.  It  can  be  shown®  that  the  auto-ambiguity 
function  of  the  input  to  each  ear  can  be  estimated  from  the  spectro- 
gram of  the  signal  (i.  e.  , the  output  of  a model  for  the  peripheral 
auditory  system).  Given  the  two  resulting  auto-ambiguity  functions 
Xuiuj(T»  and  XU2U2(T»  we  would  like  to  construct  the  magnitude- 
squared  cross -ambiguity  function  0)1^*  In  man-made  sys- 

tems, this  construction  can  be  accomplished  by  making  use  of 

9 

Stutt's  Fourier  transform  theorem: 


lXUlu2<-T*  0>|Z  = F_1 


-j]  s«1Mi»24MeiI,|>H 


Tf) 


dtdf 


where  the  operation  on  the  right-hand  side  is  a two-dimensional 
inverse  Fourier  transform  of  the  product  of  the  two  auto-ambiguity 
functions. 


Using  Stutt's  Theorem,  it  is  straightforward  to  form  a cross- 
ambiguity function  for  simultaneous  estimation  of  interaural  delay 
(angle)  and  interaural  frequency  shift  (tangential  velocity  or  evasive 
maneuver). 


r 
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The  cross-ambiguity  function  predicts  the  behavior  of  a 
system  that  uses  cross-correlation  to  measure  interaural  velocity 
difference.  Such  a prediction  can  be  checked  against  the  perfor- 
mance of  an  animal  in  an  appropriate  discrimination  experiment.  If 
performance  and  prediction  are  closely  matched,  then  one  has  support 
for  the  hypothesis  that  tangential  velocity  is  estimated  by  a cross- 
correlation process.  A similar  approach  has  been  used  by  Sayers 
and  Cherry,4  who  have  shown  that  a listener's  uncertainty  about 
sound  direction  corresponds  to  ambiguities  in  the  interaural  cross- 
correlation function.  The  method  has  also  been  used  by  Simmons  to 
compare  bat  range  resolution  performance  with  the  performance  of  a 
correlator  that  is  matched  to  the  bat's  echolocation  signal/ 

3.3.2  Standard  Deviation  of  a Cross-Range  Velocity 
Estimate 

The  following  calculations  demonstrate  the  feasibility  of 
tangential  velocity  estimation  in  animal  echolocation. 

Consider  two  transducers  that  are  d meters  apart,  and  that 
are  each  r meters  from  a target  (Figure  1). 


The  target  has  velocity  components  and  which  correspond  to 
radial  and  tangential  velocities,  respectively. 

The  radial  velocity  at  ear  #1  is 


r,  = v cos  0 + v,  sin0 
1 r t 


and  the  radial  velocity  at  ear  #2  is 


r_  = v cos 0 - v sin0, 
2 r t 


■where 


sin0  = d/Zr 


From  (1)  — (2),  we  have 


r - r_  = 2 v sin  0 
1 2 t 


Consider  a bat  such  as  Rhinolophus  and  a moth  that  is  undertaking 
evasive  action  with  tangential  velocity  v^.  We  assume  that  it  is 
sufficient  for  the  bat  to  estimate  the  direction  of  the  moth’s  motion, 
i.  e.,  the  sign  of  v^.  In  this  case,  it  is  sufficient  to  require  that 
the  uncertainty  in  measuring  v^  is  less  than  jv^j,  the  magnitude  of 
the  tangential  velocity  component.  In  other  words,  we  require 


Avt  (5) 

where  Av^  is  the  standard  deviation  of  the  velocity  measurement. 
For  a narrowband  signal  and  a white,  Gaussian  noise  background, 
an  optimum  velocity  estimator  has  standard  deviation10 
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(6) 


2 


Ap  = (SNR)  2 T 


where  T is  the  observation  time  (the  length  of  the  bat's  pulse)  and 
SNR  is  signal  energy  divided  by  noise  power  spectral  density  (signal 
to  noise  ratio).  The  quantity  A0  in  (6)  is  the  uncertainty  in  Doppler 
shift  p,  where 


= 2f  f/c 
o 


(7) 


and  where  f is  the  center  frequency  of  the  signal,  r is  range  rate, 
o 

and  c is  the  speed  of  sound.  In  our  case,  r is  either  r or  r , 

1 A 

f ~ 80  kHz,  and  c ~ 330  m/s  . 


The  uncertainty  of  is,  from  (4), 


Avt  = Affj  - r / 2 sin  6 


(8) 


where,  if  the  noise  processes  at  ear  1 and  ear  2 are  independent, 
A(r j - *2)  = Afj)2  + (Af2)2  = 42  Afj  (9) 


and,  from  (7) 


Ar  = c A*  /2fo 


Combining  (3)  — (10), 


(10) 


Avt  = rc/(<i2  SNR  TfQd). 


(11) 


Suppose  that  Rhinolophus  first  detects  its  prey  at  r = 4 m.  At  this 
range,  we  assume  that  SNR  = 3dB.  At  r = 50  cm,  SNR  is  then 

4 

greater  than  2x8  , where  the  acoustic  absorption  of  the  air  has 
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been  neglected.  It  follows  from  (11)  that,  at  r = 50  cm  and  for 
d = 2 cm, 

Avt  < (13T  x 102)"1  . (12) 

An  average  Rhinolophus  cruising  pulse  is  67  msec  long,11  but  the 

pulses  shorten  as  range  decreases.  If  T = 2.  5 msec,  Av  <0.3 

m/sec.  This  resolution  should  be  sufficient  to  determine  the 

direction  of  most  evasive  maneuvers. 

Wideband,  Doppler  tolerant  signals  are  used  by  some  FM 

^ats  ^ and  by  cetaceans.  3’  Such  signals  can  be  used  for  velocity 

estimates  if  range  rate  is  measured  with  more  than  one  echo  pulse. 

When  multiple  pulses  are  used,  it  would  seem  that  tangential 

velocity  should  be  estimated  by  measuring  rate  of  angle  change, 

i.e.,  difference  in  estimated  target  direction  for  two  different  echo 

pulses.  One  component  in  an  angle  measurement  is  the  relative 

delay  between  the  ears,  which  is  proportional  to  r - r . The 

1 2 

rate  of  angle  change  is  then  dependent  upon  measurement  of  the  time 
derivative  of  (r  - r ),  i.e.,  upon  r -r  . From  (4),  r,  - f „ can 
be  interpreted  as  a tangential  velocity  measure.  The  uncertainty 
Avt  is  thus  partially  indicative  of  the  uncertainty  in  an  angle  rate 
measurement. 

For  a train  of  wideband  pulses,  range  rate  induces  a scal- 
ing effect15  and  the  pulse  train  is  stretched  or  compressed  by  a 
factor  s,  where 

s « 1 + 2f/c  . (13) 

The  uncertainty  in  measuring  s depends  upon  the  time-bandwidth 
product  TW  of  the  pulse  train16,17 


(14) 


As  = (SNR)“2  (TW)  . 


From  (13)  - (14), 


Af  = c (SNR)  2 (ZTW) 


Using  (8)  - (9),  we  have 

_i.  -1 

Avt  = rc  (2  SNR)  2 (TW  d)  . (16) 

Comparison  of  (1 6)  and  (11)  shows  that  center  frequency  f for 
narrowband  measurement  of  Av  is  replaced  by  bandwidth  W for 
wideband  measurement  with  a pulse  train.  If  a new  pulse  is  trans- 
mitted shortly  after  a previous  pulse  is  received,  then  T « 2Nr/c, 
where  N is  the  number  of  pulses  that  are  used  for  the  velocity 
estimate,  and 


? --  -1 

Avt  = c (2  SNR)  2 (2NW  d)  . (17) 

For  the  FM  bat,  Natalus  mexicanus,  we  have  W » 115  kHz 

if  harmonics  are  used.  If  initial  detection  takes  place  at  SNR  = 3 dB 

and  at  4 meters  from  the  bat,  then  (neglecting  energy  absorption  by 

4 

the  air)  the  SNR  at  r = 50  cm  is  2 x 8 . Again  assuming  that 
d = 2 cm,  and  c = 330  m/s,  we  have,  for  r = 50  cm. 


Avt  = 1 / (6N)  m/s  . (18) 

Two  pulses  (N  = 2)  are  thus  sufficient  to  obtain  an  accurate  estimate 
of  tangential  velocity  for  assessment  of  evasive  maneuvers  at 
r = 50  cm. 

For  the  Atlantic  bottlenose  dolphin,  Tursiops  truncatus, 
we  have  W « 100  kHz,  d « 20  cm,  and  c » 1500  m/s.  Detection 


18 

distances  for  the  dolphin  are  larger  than  for  the  bat,  and  we  can 
assume  that  SNR  = 3 dB  at  least  10  m away  from  the  animal.  For 

4 

r = 1 m,  we  then  have  SNR  = 2 x 10  , and  (17)  gives 

Av^  = 1/4N  m/s  (19) 

for  r = 1 m.  Two  pulses  (N  = 2)  are  again  sufficient  for  accurate 
estimation  of  tangential  velocity  for  dolphin  echolocation. 


3.3.3  Experiments  to  Test  Sensitivity  to  Tangential 
Velocity  in  Mammals 


The  simplest  experiment  is  to  test  whether  an  illusion  of 
movement  can  be  caused  in  a human  listener  by  introducing  a time 
gated  sinusoid 


Uj(t)  = sin  ( 2 7T fQt)  rect  (t)  (20) 

into  one  ear  and  another  gated  sinusoid  of  slightly  different 
frequency 

u->(t)  = sin  [ 2 7r  ( f + 0)  t]  rect  (t)  (21) 

Ct  O 

into  the  other  ear  by  means  of  earphones.  In  (20)  - (21),  rect  (t) 
is  a rectangular  on-off  function  or  time  gate.  The  direction  of  the 
apparent  movement  should  depend  upon  the  sign  of  0,  and  the  speed 
of  movement  should  depend  upon  its  magnitude. 

A more  complicated  experiment  is  to  construct  signals  that 
have  deliberate  velocity  ambiguities,  such  as  the  sum  of  the  two 
signals  in  (20)  and  (21).  It  should  be  theoretically  possible  to 
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correlate  listener  confusion  about  apparent  motion  with  local  peaks 
in  the  cross-ambiguity  function  of  the  signals  that  are  supplied  to 
each  ear. 

For  echolocating  animals,  sensitivity  to  tangential  velocity  can 
be  tested  by  using  an  array  of  loudspeakers  that  are  placed  in  a 
semicircle  around  the  animal,  or  by  mechanically  moving  a sound 
source  or  target  on  a semicircular  trajectory.  If  the  animal  is  trained 
to  signify  the  direction  of  movement,  then  the  limits  of  its  measure- 
ment ability  can  be  tested  for  very  slow  motion,  and  resolution  can  be 
tested  by  using  two  sound  sources  that  move  at  slightly  different 
rates. 

A difficult  aspect  of  the  resolution  experiment  is  to  distin- 
guish between  a system  that  measures  rate  of  angle  change  and  a 
system  that  resolves  two  motionless  targets  at  different  angles. 

One  approach  is  to  determine  an  angle-resolution  cell  such  tha  two 
small  motionless  targets  cannot  be  resolved  if  they  are  both  within 
the  cell,  and  to  perform  all  moving  target  experiments  with  objects 
that  are  not  separated  by  more  than  a single  angle  resolution  cell. 

For  targets  that  move  at  a constant  rate,  the  observation  time  must 
be  limited  to  the  time  when  the  targets  are  close  together. 

3.3.4  Conclusion 

In  order  to  compensate  for  evasive  maneuvers  at  close  range, 
a sonar  should  be  capable  of  determining  the  direction  of  cross- 
range movement  when  the  sonar  is  close  to  the  target.  This  capa- 
bility exists  in  dolphin  echolocation  if  two  succeeding  pulses  can 
be  coherently  processed.  A few  more  pulses  may  be  necessary  if 
suboptimal  processing  is  used. 

For  some  man-made  sonar  systems,  a cross-range  velocity 
estimate  is  relevant  not  only  to  measurement  of  evasive  maneuvers 
but  also  to  the  use  of  target  motion  as  a clutter  rejection  method. 
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Measurement  of  cross-range  velocity  is  also  important  for  track- 
ing and  beam-forming.  The  best  boresight  position  for  a future 
transmission  can  be  predicted  if  a cross-range  velocity  estimate 
is  available. 

In  the  above  applications,  one  would  want  to  measure  cross- 
range velocity  at  a long  distance  from  the  target.  This  measurement 

1/2 

can  be  accomplished  by  increasing  SNR  Tf  d in  (]l),  or  by  increas- 

1/2  ° 

ing  SNR  TWd  in  (16). 

1/2 

To  increase  SNR  , one  can  use  pulse  compression  tech- 
niques with  long  signal  durations.  To  increase  d,  one  can  gather 
data  from  two  ends  of  a large  towed  array  or  (better  still)  from 
helicopter-dipped  transducers  that  communicate  with  a central 
processor.  For  long  target  ranges,  the  signal  duration  T can  be 
significantly  increased.  An  increase  in  bandwidth  W may  be  prefer- 
able to  an  increase  in  f^,  since  attenuation  increases  with  frequency. 

The  same  bandwidth-array  size  trade-off  that  was  found  in 
Section  3.  2 for  target  position  measurement  is  again  found  in  (16). 
Both  angle  and  angle  rate  measurements  can  be  improved  by  in- 
creasing array  size  or  by  increasing  the  bandwidth  of  the  sonar 
system. 
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3.4 


A Test  of  the  Energy  Spectrum  Analysis  (Johnson-Titlebaum) 
Hypothesis 

3.4.1  Introduction  to  Section  3.4 


R.  A.  Johnson  and  E.  L.  Titlebaum  (1976)  have  proposed  a 
model  for  echolocation  processing  that  is  related  to  the  perception 
of  time  separation  pitch.  The  proposed  model  forms  the  magnitude- 
squared  Fourier  transform  of  the  transmitted  pulse  and  an  echo, 
i.  e.  , transmitted  and  echo  pulses  are  both  included  in  the  integra- 
tion interval  for  Fourier  analysis.  For  an  echo  with  delay  A and 


attenuation  a,  the  resulting  energy  density  spectrum  is 


|R(f)|2  = |U(f)|2  [1  + a2  + 2 acos  (27rAf)] 


where  U(f)  is  the  Fourier  transform  of  the  transmitted  signal.  The 
presence  of  a target  introduces  a frequency  domain  ripple  with  period 
1 /A  Hz.  Detection  of  a target  and  estimation  of  its  range  depend 

upon  an  analysis  of  this  ripple. 

In  his  thesis,  R.  A.  Johnson  (1972)  has  suggested  that  rippled 

noise  could  be  very  detrimental  to  the  performance  of  an  energy 
spectrum  analyzer.  In  this  section,  we  use  locally  optimum  detection 
ideas  to  provide  further  insight  into  the  effect  of  rippled  noise. 


3.4.2  Locally  Optimum  Detection  Using  the  Energy 
Spectrum  as  Data 

Our  objective  is  to  obtain  a statistical  test  that  can  distinguish 
between  the  hypothesis  Hj  that  a target  is  present  in  additive  Gaussian 
noise  and  the  hypothesis  that  only  Gaussian  noise  is  present.  Let 
P (f)  be  the  noise  power  spectral  density.  The  expected  form  of  the 
energy  spectrum  under  the  two  hypotheses  is  then 


Ho:  E 


{ I R(f)  + N(f)  I2}  = |u(f)|2+PNJ 


H . : e|  | R(f)  + N(f)  ] = ]u(f)  |2[l+a2  +2acos  (2ffAf)]  + P (f). 

1A  * 

where  H is  the  hypo^.sis  that  a target  is  present  at  range  ca/2, 
and  N(f)  is  the  Fourier  transform  of  a sample  function  of  the  noise 
process. 

One  can  obtain  a set  of  uncorrelated  samples  of  the  energy 
spectrum  by  using  samples  that  are  (4T)  *Hz  apart,  where  T is  the 
integration  time  of  the  Fourier  analyzer.  Let  z^  denote  the  sample 
at  frequency  f . The  probability  densities  of  z.  under  the  two  hypo- 

j J 

theses  are  then 


z.|Hq)  = Cr."2exp|[-z.-  | U(fj)l2]  / 


• I0[2x/du(f.)  I a.2] 


p(z.\  Hj)  = a.  2 exp  |{”zj“  lU(fj>  I2!1  + ^ + 2acos  ( 2 TT  A f / a2j 
1^2  4|u(f.)  I 2[  1 +a2+  2a  cos(2J7Af.)]  /a2} 


where  <7.  = P_T(f. )• 

J N j 

We  now  define  a signal-to-noise  ratio  parameter  6,  where 

e = a/a ^ 

Replacing  IQ(-)  by  its  power  series  representation,  we  have 
p(z.|H0)  = <7j2  exp|f-z.  - |U(fj  )|2]/a2} 

[1+  \ U(f.)  1 2z./a4  + ...] 

J J J 


i 
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p(z  |H,  ) = O . exp  <-cr.  z.-|U(f.)|[cr.  + 0 + 2ff  . ecos(2ffAf 

r 3 \ J J J J J 

• jl  + 0~2  z.  |U(f.)  |2[a.2  + e2  + 2a.10cos(2ff  M.)]  + ...  | 
J J J J J J 

The  log  likelihood  ratio  for  M independent  samples  is 

M M 

In  A = ln[  17  p(z.  H )/ IT  p(z.  H )]. 

j=l  J1  j=l  J u 

The  locally  optimum  test  for  deciding  between  H^and  HQ  (Capon, 
1961;  Middleton,  1966)  is  to  compare 
(d/d9)  In  I 


with  a threshold,  where 

M 

(d/de)  In  aJ  = £[-2a.|u(f.)|  cos(2rrAf.)] 

>0  J=1  J J J 

M [2a73  z.|U(f.)|2cos(27TAf.)] 

+ E — 1 — L-r1 j — ^ ' 

j=1  i + cr  z |u(f  )I 

The  part  of  (5)  that  is  independent  of  the  data  z.  can  be  incorporated 

into  a threshold  level.  The  locally  optimum  test  is  then  to  compare 
M 

l(z)  = E a.  z.  cos(2ffM.)  (6) 

j=l  J J J 

with  a threshold,  where 

a = 2a  3 lu(fjl2/[l+a'4z.|u(f;)|2] 


The  total  number  M of  frequency  domain  samples  increases 
with  signal  bandwidth.  The  likelihood  ratio  thus  involves  two  non- 
central chi-square  distributions  with  degrees  of  freedom  that  depend 
upon  signal  bandwidth.  The  locally  optimum  detector  in  (6)  (and  the 
maximum  likelihood  estimator  for  the  delay  parameter  A)  will  in- 
volve correlation  of  the  sampled  energy  spectrum  with  a vector 
that  exhibits  a hypothetical  ripple,  i.  e. , Fourier  analysis  of  the 


energy  spectrum. 


If  P (f)  in  (2)  has  a ripple  similar  to  that  in  (1),  then  similar 
N 


energy  density  spectra  will  be  observed  under  both  hypotheses  and 
it  will  be  difficult  to  distinguish  between  then. 


3.4.3  Effect  of  Rippled  Noise  Upon  a Matched  Filter 


A matched  filter  has  an  expected  output  noise  power  P , 


where 


00 

/pN(f)|U(f)  I2  df 


under  both  hypotheses.  A ripple  in  P (f)  will  not  greatly  affect 


P if  | U(f)  | is  smooth  and  broadband  (as  is  the  case  for 
MF 


many  dolphin  clicks  and  FM  bat  pulses).  Specifically,  a periodic 


ripple  in  P^(f)  with  period  1/  A Hz  will  have  very  small  effect  upon 


the  integral  in  (7)  if 


1/A  <<  B 


where  B is  the  signal  bandwidth.  An  adaptable  whiten-and-match 
u 


filter  will,  in  fact,  take  advantage  of  the  nulls  in  a rippled  noise 


power  spectral  density  P^.(f)  in  order  to  improve  detection  perfor- 


mance . 


......  . . ...  ......  . . . 
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3.4.4  Experimental  Procedure  for  Testing  the  Johnson-Titlebaum 
Model  Against  the  Matched  Filter  Model 

Assuming  that  the  Johnson-Titlebaum  model  is  implemented 
by  an  echolocating  animal,  broadband  masking  noise  will  be  espe- 
cially effective  if  the  noise  has  power  spectral  density 

pN(f)  = <N0/2)  [1  + a2  + 2ttcos  (27rAf)]  (9) 

where  A is  the  delay  that  is  associated  with  the  target.  Masking 
noise  with  the  power  spectral  density  in  (9)  should  have  a much 
greater  effect  upon  the  energy  spectral  derfsity  model  than  white 
noise  with  the  same  average  power.  The  appropriate  noise  process 
can  be  generated  by  passing  white  noise  through  an  attenuator  and 
a delay  line,  and  by  adding  the  attenuated,  delayed  noise  to  the 
original  noise  process. 

The  noise  can  be  transmitted  toward  the  animal  from  behind 
a small  target  which  the  animal  attempts  to  detect.  Detection  per- 
formance should  be  maximally  degraded  when  A in  (9)  is  equal  to  the 
delay  between  transmitted  pulse  and  echo,  and  there  should  be  less 
degradation  for  other  values  of  A. 

If  the  matched  filter  model  is  correct,  then  performance 
degradation  should  be  independent  of  the  parameter  A in  (9),  pro- 
vided that  condition  (8)  holds. 

The  same  masking  noise  can  be  used  for  human  subjects  who 
are  attempting  to  estimate  time  separation  pitch  or  to  discriminate 
a single  pulse  from  a pulse  pair,  where  the  second  pulse  has  smaller 
amplitude  than  the  first.  These  tests  can  be  used  to  verify  our  pre- 
diction that  a rippled  noise  spectrum  is  especially  disturbing  to  a 
detection  and  range  estimation  process  that  uses  energy  spectral 
analysis . 
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A reasonable  experimental  procedure  has  two  steps: 

1.  Given  two  broadband  pulses  u(t)  + au(t  - Aj ),  find  the 
threshold  value  of  a in  human  subjects  such  that  the  second 
pulse  is  detected  (say)  75%  of  the  time  in  a white  noise  back- 
ground. Call  this  value  a . Also  find  the  threshold  value 
of  a for  detection  of  the  second  pulse  in  a rippled  noise 
background  with  (f ) as  in  (9),  such  that  P^(f)  has  the 
same  average  power  as  the  white  noise.  Call  this  value 
Evaluate  the  effect  of  rippled  noise  upon  human  estimation 
of  time  separation  pitch  when  a = 1 . 

If  a > oc  and/or  if  the  rippled  noise  in  (9)  maxi- 
r w 

mally  degrades  the  estimation  of  time  separation  pitch  when 
A = Aj,  proceed  to  the  following  experiment. 

2.  Evaluate  the  effect  of  masking  noise  with  rippled  power  spec- 
tral density  as  in  (9)  upon  target  detection  and  range  estima- 
tion performance  of  dolphins  and  bats.  Use  white  noise  with 
the  same  average  power  and  a sequence  of  different  A-values 
in  (9)  for  comparison. 

3.4.5  Conclusion 

Both  intuition  and  a locally  optimum  detector  formulation 
indicate  that,  to  detect  a target  from  an  energy  spectrum  that  in- 
cludes the  transmitted  signal  as  well  as  hypothetical  echoes,  one 
should  perform  a Fourier  analysis  of  the  energy  spectrum.  Such 
an  analysis  is  susceptible  to  rippled  noise,  i.  e.  , a noise  power 
spectrum  that  is  modulated  with  a periodic  function.  If  the  period 
of  this  ripple  is  small  relative  to  signal  bandwidth,  then  the  ripple 
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will  not  greatly  affect  the  performance  of  a matched  filter.  This 
analysis  leads  to  a straightforward  behavioral  test  that  can  be  used 
to  decide  whether  matched  filtering  or  energy  spectrum  analysis  is 
a better  model  for  animal  echolocation . 
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